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METHOD AND APPARATUS FOR GENERATING AND DISTRIBUTING 
SATELLITE TRACKING INFORMATION IN A COMPACT FORMAT 

BACKGROUND OF THE INVENTION 
Field of the Invention 

[0001] The present invention generally relates to generating satellite tracking 
information for earth orbiting satellites. More specifically, the invention relates to a 
method and apparatus for generating and distributing satellite tracking information in a 
first format (e.g., a compact ephemeris model) through a network or communications 
link, then representing the satellite tracking information in a second format (e.g., a 
standard ephemeris model) at a receiver. 

Description of the Related Art 

[ooog A positioning receiver for the Global Positioning System (GPS) uses 
measurements from several satellites to compute a position. The process of acquiring 
the GPS radio signal is enhanced in speed and sensitivity if the GPS receiver has prior 
access to a model of the satellite orbit and clock. This model is broadcast by the GPS 
satellites and is known as ephemeris or ephemeris information. Each satellite 
broadcasts its own ephemeris once every 30 seconds. Once the BPS radio signal has 
been acquired, the process of computing position requires the use of the ephemeris 
information. 

pop's] The broadcast ephemeris information is encoded in a 900 bit message within 
the GPS satellite signal. It is transmitted at a rate of 50 bits per second, taking 18 
seconds in all for a complete ephemeris transmission. The broadcast ephemeris 
information is typically valid for 2 to 4 hours into the future (from the time of broadcast). 

Before the end of the period of validity the GPS receiver must obtain a fresh 
broadcast ephemeris to continue operating correctly and produce an accurate position. 
It is always slow (no faster than 18 seconds), frequently difficult, and sometimes 
impossible (in environments with very low signal strengths), for a GPS receiver to 
download an ephemeris from a satellite. For these reasons it has long been known 
that it is advantageous to send the ephemeris to a GPS receiver by some other means 
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in lieu of awaiting the transmission from the satellite. US Patent 4,445,1 18, issued April 
24, 1984, describes a technique that collects satellite orbit information at a GPS 
reference station, and transmits the information to the remote GPS receiver via a 
wireless transmission, This technique of providing the ephemeris, or equivalent data, 
to a GPS receiver has become known as "Assisted-GPS". Since the source of 
ephemeris in Assisted-GPS is the satellite signal, the ephemeris information remains 
valid for only a few hours. As such, the remote GPS receiver must periodically connect 
to a source of ephemeris information whether that information is received directly from 
the satellite or from a wireless transmission. Without such a periodic update, the 
remote GPS receiver will not accurately determine position. 

[0004] Furthermore, the Assisted-GPS systems typically retransmit the entire 
ephemeris message to the remote receiver. In many instances, bandwidth or packet 
size for the transmission of this message is not readily available. 

10005] Therefore, there is a need for a method and apparatus for providing satellite 
trajectory and clock information to a remote receiver in a compact form. 

■SUMMARY OF THE INVENTION 

[0006] The present invention is a method and apparatus for generating satellite 
tracking data (STD), then transmitting the data to a remote receiver in a compact form. 

The STD is derived by receiving at one or more satellite tracking stations the signals 
from at least one satellite and determining satellite tracking information (STI) through 
signal processing or by extracting the ephemeris message from the received signals. 
STI contains present satellite orbit trajectory data and satellite clock information. 

[00 07] The STD is reformatted into a compact format and provided to a remote 
satellite signal receiver via a network or communications system. The receiver 
converts the compact format into a standard format and uses the STD to compute the 
position of the receiver. The satellite system may include the global positioning system 
(GPS), GLONASS, GALILEO, or other satellite systems that may use STD to enhance 
the performance of the receiver. 
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BRIEF DESCRIPTION OF THE DRAWINGS 
[ooos] So that the manner in which the above recited features of the present 
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invention are attained and can be understood in detail, a more particular description of 
the invention, briefly summarized above, may be had by reference to the embodiments 
thereof which are illustrated in the appended drawings. 

[0009] It is to be noted, however, that the appended drawings illustrate only typical 
embodiments of this invention and are therefore not to be considered limiting of its 
scope, for the invention may admit to other equally effective embodiments. 

{0010] Figure 1 depicts a system for creating and distributing satellite tracking data 
(STD) to remote GPS receivers; 

rooii] Figure 2 depicts a flow diagram of a method for forming the STD from the 
satellite measurements made at satellite tracking stations; 

[0012] Figure 3 depicts a flow diagram of a method for forming a compact orbit 
model in accordance with the present invention; and 

[0013] Figure 4 depicts an example of compacting the orbit model, where two orbit 
model terms ere compacted into a single term. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 

[0014] Figure 1 depicts a block diagram of a system 100 for creating and distributing 
satellite tracking data (STD). The satellite system may include the global positioning 
system (GPS), GLGNASS, GALILEO, or other satellite systems that may use STD to 
enhance the performance of the receiver, The following disclosure uses GPS as an 
illustrative system within which the invention operates. From the following disclosure, 
those skilled in the art will be able to practice the invention in conjunction with other 
satellite based positioning systems. 

[0015] A network of GPS tracking stations 102 is used to collect measurement data . 
from the GPS satellites 104. Such a network is described in detail in US patent 
application serial number 09/615,105, filed July 13, 2000 and incorporated herein by 
reference. The network could comprise several tracking stations that collect satellite 
tracking information (STI) from all the satellites in the constellation, or a few tracking 
stations, or a single tracking station that only collects STI for a particular region of the 
world. An STD collection and computation server 1 06 collects and processes the 
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measurement data (this measurement data is referred to herein as .satellite tracking 
information (STI)that includes at least one of; code phase measurements, carrier 
phase measurements, Peppier measurements, or ephemeris data). The ephemeris 
data may be the decoded ephemeris message extracted from the GPS signal itself. 

[0016] The server may create long term STD in accordance with the teachings of 
U.S. patent application serial number 09/875,809, filed June 6, 2001, (attorney docket 
GLBL/016) or standard ephemeris message data in accordance with the teachings of 
U.S. patent number 5,365,450, issued November 15, 1994, both of which are 
incorporated herein by reference. The server 106 may produce one or more of the 
following: 1) accurate satellite tracking data (STD) (e.g., a trajectory of each satellite 
and/or a clock offset measurement) during the data collection period, 2) a prediction of 
the future STD of each satellite, and 3) models that match the future STD of each 
satellite. 

rooiT] The server 106 comprises a central processing unit (CPU) 118, support 
circuits 122, and memory 120. The CPU 118 may be any one of the many CPUs 
available on the market to perform general computing. Alternatively, the CPU may be 
a specific purpose processor such as ah application specific integrated circuit (AS! C) 
that is designed to process satellite tracking information. The support circuits 122are_ 
well known circuits such as clock circuits, cache, power supplies and the like, Trie 8 ** 
memory 120 may be read only memory, random access memory, disk drive storage, 
removable storage or any combination thereof. The memory 120 stores executable 
software, e.g., STD software 124, thai, when executed by the CPU 118, causes the 
system 100 to operate in accordance with the present invention. 

[0018] The set of satellite trajectory and clock data produced by the STD software 
124. The STD is stored in an STD database 108. A distribution server 110 accesses 
the database 108 to gather the most recent set of STD, formats the data using the 
formatting software 111, and distributes the formatted data to GPS devices 1 1 2 that 
require satellite orbit information. The software 11 1 produces a compact format, e.g., £ 
compact ephemeris model, in accordance with the present invention. 

[ DD19] The distribution process may be implemented using some form of wireless 
communications system 114, or over the Internet 116, or a combination of both, or by 
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some other means of communication. Once Ihe GPS devices 112 have received, the 
compact epherneris model, they expand the model to a format that is conventional for 
receiver. The compact epherneris mode! distributed to the GPS devices may be in a 
similar format as the broadcast epherneris or may be some other model format that is 
defined by the GPS device. Herein this orbit data is generally referred to as a satellite • 
tracking model (STM). The loading of the STM into the GPS receiver can be 
accomplished in many ways. Using the cradle for a personal digital assistant (PDA), 
direct connection to a network, or a wireless technology, such as Bluetooth or a cellular 
network, are a few examples of how tie satellite data can toe transferred to the 
receiver. The transmission is generally accomplished by broadcasting a compact 
model of the STD (or a compact model representing a portion of the STD) without 
knowledge of the specific location of the GPS receiver. As such, the distribution .server 
does not require the GPS receiver to send any information through the network to ike 
distribution server. 

p)2p] Since GPS is a ranging system in land of itself, the data transmitted toy the 
GPS satellites can be used to determine the range, range-rate and clock offsets to the 
GPS satellites from a set of tracking stations. This set of observations generated bythe 
tracking stations 102 is used in the orbit determination process, and in the estimation 
of the satellite clock characteristics. The set of monitoring stations 102 could be a 
single station, a public network such as the Continuously Operating Reference System 
(CORS), or a privately owned and/or operated network. 

[0021] FIG. 2 depicts a flow diagram of the process 20G of the present invention. 
The process 200 begins at step 202, wherein the satellite measurements are collected 
at the tracking stations. At step 204, the satellite trajectory data (STD) is computed or 
extracted from the satellite signals. The STD is then stored at step 206 in the STD 
database. At step 208, the database is accessed and the formatting software is 
executed to convert the formatting of the accessed STD, The formatted STD is output 
as the compact model at step 210. 

[0022] One embodiment of the invention formats the STD as a subset of the 
standard epherneris parameters defined in ICD-GPS-200c. Fitting the STD to the 
desired compact orbit model can be accomplished in a number of mathematical 
methods. The preferred embodiment is a least-squares fit of the orbit model 
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parameters to the trajectory data. Other methods, such as Kalman filters or other 
estimators can also be used to obtain the orbit model parameters that best fit the 
trajectory data. These techniques of fitting data to orbit models are well known to 
people skilled in the art of orbit determination and orbit modeling. 

[0023] The least squares technique provides an optimal fit of the trajectory data to 
the model trajectory formed from the compact orbit model parameters. Figure 3 
depicts a flow diagram of a method of generating an orbit model using a least squares 
estimation technique. 

'■ [0024] At step 302, the STD for the desired time interval is extracted from the STD 
database. The orbit model parameters are initialized to the orbit mode! values obtained 
by a similar process for the previous interval. This guarantees that the initial orbit 
model parameters are a good tit at least for the beginning of the desired time interval. 
The rest of the process 300 will ensure that the parameters are adjusted . so that they 
become a good fit for the entire time Interval, 

[0025] In the preferred embodiment there are 15 orbital parameters to be adjusted: 

• Square root of semi-major axis (meters A 1/2) 

• Eccentricity (dimensionless) 

• Amplitude of sine harmonic correction term to the orbit radius (meters) 

. Amplitude of cosine harmonic correction term to the orbit radius (meters) 
. Mean motion difference from computed value (radia ns/sec) 



« Amplitude of cosine harmonic correction term to the argument of latitude 
(radians) 

• Amplitude of sine harmonic correction term to the argument of latitude 
(radians) 

o Amplitude of cosine harmonic correction term to the angle of inclination 
(radians) 

s Amplitude of sine harmonic correction term to the angle of inclination 
(radians) 

. Longitude of ascending node of orbit plane at weekly epoch (radians) 
m Inclination angle at reference time (radians) 

6 
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. Rate of inclination angle (radians/sec) 

o Argument of perigee (radians) 

» Rate of right ascension (radians/sec) 

P027] At step 303, some of the terms in the 15 term set are set to zero. The terms 
that are selected are the 6 harmonic terms such that there are 9 remaining parameters. 
This approach is particularly useful when bandwidth and/or packet size is limited in the 
communication link that Will be used to convey the orbit model to the satellite Signal 
receiver, e.g., the remote GPS receiver. The subset Of 9 parameters, by setting all 
harmonic terms in the model to zero, is: 

. Square root of semi-major axis (meters A 1/2) 
. Eccentricity (dimensionless) 

. Mean motion difference from computed value (radians/sec) 
. Mean anomaly at reference time (radians) 

. Longitude of ascending node of orbit plane at weekly epoch -(radians) 

* Inclination angle at reference time (radians) 

. Rate of inclination angle (radians/sec) 

o Argument of perigee (radians) 

. Rate of right ascension (radians/sec) 
:f0028] the receiver can then reconstruct a standard epheme'ris model by setting the 
"missing" harmonic terms to zero. In essence, the receiver reformats the -STD for 
processing by the receiver circuits. 

po2Si] As an example of the method of generating the compact model, consider Figure 
4, which shows, for simplicity, just two terms of an orbit 400; an orbital radius (A), and a 
radial harmonic term (r). For this simple example, these two terms form the non- 
compact model, wherein the orbit is described by a circle of radius (A) plus a harmonic 
perturbation (r). To produce a more compact model that fits the actual orbit over an 
interval 402, the method of the invention removes the harmonic term (i.e., sets the term 
(r) to zero) and increases the orbital radius (A) to a larger value (A,). The compact 
model is an orbit described by a circle with radius A-,. If an application requires a non- 
compact orbit model, then the compact model (Ai) can be represented as a non- 
compact model by specifying a harmonic term (r, ) equal to zero. This compact model 
7 
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will fit the original orbit, over an interval 402, with a small error. 

[0030] to the preferred embodiment, 6 harmonic terms are removed from trie 15- 
parameter model, and the other 9 terms are adjusted by process 300 that is analogous 
to the example 400 to provide a compact model that is accurate over a pre-defined 
interval. By adjusting the 9 remaining terms of an orbit model, while "zeroing" 6 
harmonic terms, the compact model can be made accurate over a period of time such 
that a GPS receiver that relies on a compact mode! to compute position would 
compute a location that is no more than 2 meters less accurate than if the receiver 
used a full orbit model to compute position. 

[0031] There are many alternative embodiments that will be readily apparent to those 
skilled in the art, such as removing more or fewer terms before adjusting the remaining 

terms, setting removed terms to some value other than zero, and defining hew terms ' 

that model the orbit. 

[0032] Returning to FIG. 3, at step 304, the orbit model is used to predict what the 
trajectory would be, the predicted data is denoted the "Model Trajectory Data" (MTD). 
If the model were perfect, the MTD would coincide exactly with the STD. At step 306, 
the MTD and STD are compared to see how closely the orbit model fits the orbit data. 
In the preferred embodiment, the comparison step 306 is performed by summing the 
squares of the differences between each trajectory point in the STD and the c 
corresponding point in the MTD, and comparing the resulting sum to a threshold. If the 
fit is "good", the model parameters are deemed "good" and the process stops at step 
310. If the fit is not good then the model parameters are adjusted at step 308. There 
are many techniques well known in the art for adjusting model parameters to fit data. 
Steps 304, 306 and 308 are repeated until the model parameters are found that fit the 
STD well. 

£0033] There are a large number of alternative embodiments to reduce the size of the 
data, i.e., compacting the STD, while still providing a model that fits the STD, including: 

* Removing parameters from the model, and replacing them with a 

constant, such as zero - as done above - or some other predetermined =t 
value, which is either stored in the Remote GPS Receiver, or occasionally 
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sent to the receiver. The predetermined value may be determined by a 
GPS almanac stored at both the receiver, and the distribution server. 
» The resolution of the parameters may be restricted in the process 300, 
this too reduces the amount of data that must be sent to the mobile GPS 
receiver. 

- Parameters, which are similar among two or more satellites, may be 
represented as a master value plus a delta, where the delta requires 
fewer bits to encode; an example of this is the parameter Eccentricity, 
which changes very little among different GPS satellites. 

[0034] Some of these approaches reduce the ability of the model to fit the data over a 
period of time (e.g., four hours). In this case, the fit interval may be reduced {e.g. to 
two hours) to compensate. The accuracy of fit of the model can be traded off against 
the period of time over which the model is valid. While the foregoing is directed to the 
preferred embodiment of the present invention, other and further embodiments of the 
invention may be devised without departing from the basic scope thereof/ and the 
scope thereof is determined by the claims that follow. 
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Claims: 

1 . A method of creating and distributing compact satellite orbit models comprising: 
receiving satellite signals from at least one satellite and at least one receiving 

station- 
extracting at least a portion of the satellite tracking data from said satellite 

signal, representing said data in a first format; 

transmitting the formatted data to a remote receiver; and 

at the remote receiver, representing said formatted data in a second format 

supported by the remote receiver. 

2. The method of claim 1 wherein said satellite tracking data comprises at least 
one of a satellite orbit model or a satellite clock model. 

3. The method of claim 1 wherein said satellite tracking data comprises at least 
one of: data representative of a satellite orbit model, orbit model, data representative 
of a satellite clock model. . 

A. The method of Claim I wherein said second format is a format that is prescribed 



5. The method of claim 1 wherein said second format comprises parameters 
defined in ICD-GPS-200. 

6. The method of claim 1 wherein said first format comprises a first set of terms 
that define a first form of orbit model, and said second format comprises a second set 
of terms that define a second form of orbit model, and said second set of terms is 
larger than said first set of terms. 

7. The method of claim 1 wherein said first format comprises a first set of terms 
that define a first form of orbit model, and said second format comprises a second set 
of terms that define a second form of orbit model, where said second set of terms 
contains said first set of terms. 
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B. The method of claim 1 wherein said first format comprises a first set Of terms 
that define a first form of orbit model and said second format comprises a second set 
of terms that define a second form of orbit model, and at least one term of said first set 
is defined as a number with lower resolution than the corresponding term in second 
set. 

9. The method or claim 1 where said data in a first format requires fewer bits to 
encode it than said data in a second format. 

10. The method of claim 1 wherein said remote receiver is a GPS receiver. 

11. The method of claim 1 wherein said remote receiver is a satellite positioning 
system receiver. 

12. The method of claim i wherein said second format is a standard format for 
providing satellite models to a global positioning system receiver. 

13. The method of claim 1 wherein said transmitting is performed using a wireless 
communications link. 

14. The method of claim 13 wherein said transmitting is done without requiring from 
said remote receiver data representative of said remote receiver's position. 

15. The method of claim 1 where an accuracy of the data in said first format is 
increased by decreasing a time interval represented by said formatted data, 

16. A method of creating and distributing a compact orbit model comprising: 
receiving satellite signals from at least one receiving station; 

extracting at least a portion of the satellite tracking data from the satellite signal, 
where said portion comprises a first number of orbit terms that define a first orbit 
model; 

formatting said portion to form formatted data having a second number of orbit 
terms that define a second orbit model, where said first number is greater than said 

11 ^ 
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second number; 

transmitting the formatted data to a remote receiver; and 

at the remote receiver, expanding the formatted data to have the first number df. 

terms. 

17. The method of claim 16 wherein said satellite tracking data comprises at least one 
of a satellite orbit model or a satellite clock model. 

18. The method of claim 1 6 wherein said satellite tracking data comprises at least 
one of: data representative of a satellite orbit model, orbit model, data representative - 
of a satellite clock model. 

19. The method of claim 16 wherein said second orbit model is a compact orbit model. 

20. The method of claim 16 wherein said expanding step results in a data format 
that is prescribed by said remote receiver. 

22. The method of claim 20 wherein the data format having the first number of 
terms Comprises parameters defined in lCD-GPS-200. 

23. The method of claim 16 where said terms in said second orbit model require 
fewer bits to encode it than said terms in a first orbit model. 

24. The method of claim 16 wherein said remote receiver is a GPS receiver. 

25. The method of claim 16 wherein said remote receiver is a satellite positioning 
system receiver. 

26. The method of claim 16 where an accuracy of the data in said second orbit 
model is Increased by decreasing a time interval represented by said formatted data \ 
defining said second orbit model. 



27. A method of creating a compact orbit model comprising: 
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providing a satellite orbit model having a fi rst set of terras; 

adjusting the first set of terms to prod uce a compact orbit model having a 
second set of terms, where a number of terms in the first set of terms is greater than a 
number of terms in said second set of terms. 

28. The method of claim 27 wherein said adjusting step further comprises: 

zeroing a plurality of terms in said first set of terms. 

29. The method of claim 28 wherein said adjusting step further comprises: 

adjusting a plurality of non-zero terms in the second set of terms in response to 
the effects of zeroing terms in the first set of terms. 

30. Apparatus for creating and distributing compact satellite orbit models comprising: 

at least one satellite signal receiver for receiving satellite signals from at least 
one satellite; 

means for extracting at least a portion of the satellite tracking data from said 

satellite signal, representing said data in a first format; 

a transmitter for transmitting the formatted data to a remote receiver; and - 
at the remote receiver, representing said formatted data in a second format 

supported by the remote receiver. 

31. The apparatus of claim 30 wherein said second format is a format that is 
prescribed by said remote receiver. 

32. The apparatus of claim 30 wherein said second format comprises parameters 
defined in ICD-GPS-200. 

33. The apparatus of claim 30 wherein said first format comprises a first set of terms 
that define a first form of orbit .model, and said second format comprises a second set 
of terms that define a second form of orbit model, and said second set of terms is 
larger than said first set of terms. 

34. The apparatus Of claim 30 wherein said transmitter is a wireless 
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